Submitted to ApJ Letters 

Preprint typeset using LAT^X style emulateapj v. 16/07/00 



(N 
O 
O 
(N 

X> 

<D 

PL. 
OS 



(N 
> 
(N 

m 

(N 
On 
O 



^ ■ 

Or 
i 1 

o ■ 



THE AGE OF THE UNIVERSE AND THE COSMOLOGICAL CONSTANT DETERMINED FROM 
COSMIC MICROWAVE BACKGROUND ANISOTROPY MEASUREMENTS 

Lloyd Knox 

Department of Physics, University of California, Davis, CA 95616, USA, email: lknox@ucdavis.edu 

Nelson Christensen 

Physics and Astronomy, Carleton College, Northfield, MN, 55057, USA email: nchriste@carleton.edu 

AND 

CONSTANTINOS SKORDIS 

Department of Physics, University of California, Davis, CA 95616, USA, email: skordis@bubba.ucdavis.edu 

Submitted to ApJ Letters 

ABSTRACT 

If f2tot = 1 and structure formed from adiabatic initial conditions then the age of the Universe, as 
constrained by measurements of the cosmic microwave background (CMB), is to = 14.0 ± 0.5 Gyr. The 
uncertainty is surprisingly small given that CMB data alone constrain neither h nor Q\ significantly. It 
is due to the tight (and accidental) correlation, in these models, of the age with the angle subtended 
by the sound horizon on the last-scattering surface and thus with the well-determined acoustic peak 
locations. If we assume either the HST Key Project result h = 0.72 ± .08 or simply that h > 0.55, 
we find J1a > 0.4 at 95% confidence — another argument for dark energy, independent of supernovae 
observations. Our analysis is greatly simplified by the Monte Carlo Markov chain approach to Bayesian 
inference combined with a fast method for calculating angular power spectra. 

Subject headings: cosmology: theory - cosmology: observation - methods: data analysis - methods: 
statistical - cosmology: cosmic microwave background, cosmological parameters, 
distance scale 



1. INTRODUCTION 

Determining the expansion age of the Universe has been 
a major goal of cosmology ever since Hubble discovered 
the expansion. Compatibility with determinations of stel- 
lar ages is an important consistency check of cosmological 
models. Traditional methods of determining the expan- 
sion age rely on Hubble constant measurements which are 
either highly imprecise, or have error budgets dominated 
by systematics. In addition one must determine f^A (or 
more generally the mean density of the various compo- 
nents) since it affects how the expansion rate has changed 
over time. In this Letter we present highly precise age de- 
terminations from CMB data which completely bypass the 
need for independent determinations of H$ and Q,\. 

There are a handful of cosmological parameters which 
can be determined from measurements of the CMB an- 
gular power spectrum to percent level accuracy, such 
as u>b, u> m , and f2tot (where a;, = flih 2 and Hq = 
lOO/ikmsec^Mpc -1 ) (e.g. Eisenstein et al. 1999). Other 
parameters can not be well-determined and require the 
addition of complementary observations. For example £l\ 
is poorly determined by the CMB alone(e.g. Efstathiou & 
Bond 1999) but well-determined when supernovae obser- 
vations are included (e.g. Netterfield et al. 2001). 

It has been pointed out (Ferreras et al. 2001) and 
demonstrated (Netterfield et al. 2001) that the CMB can 
be used to place tight constraints on the age of the Uni- 
verse. This is due to the high degree of correlation be- 
tween the angle subtended by the sound-horizon on the 
last-scattering surface, 9 S , and age in flat adiabatic mod- 
els, also noticed by Hu et al. (2001) who used it to place 
an upper bound on the age. Here we extend the previous 
work by including additional data, by taking the flatness 
assumption seriously and by use of a new analysis tech- 



nique which has advantages as described below. We also 
demonstrate the accidental nature of the age-sound hori- 
zon correlation by showing that its tightness depends on 
where we are in the h parameter space. We are fortu- 
nate that the correlation is tightest near the "concordance" 
values of h = 0.72 and = 0.65. 

Our age determination is model-dependent and the 
model (adiabatic CDM) has many parameters. We take 
them to be the amplitude and power law spectral index 
of the primordial matter power spectrum, A and n, the 
baryon and dark matter densities, tVb and u>d, the cosmo- 
logical constant divided by the critical density, and 
the redshift of re-ionization of the intergalactic medium, 
Zri. The Hubble constant and age are derived parameters, 
given in terms of the others by h 2 — (uji, + — 



and t = 6.52 Gyr ln[(l + VOvVVl - n A }/VTlJ?. 

We do not consider models with f2 t ot 7^ 1 or dark energy 
models other than the limiting case of a cosmological con- 
stant with w = Pf p = — 1. The first wc justify on grounds 
of simplicity: CMB observations indicate the mean curva- 
ture is close to zero and generally agree well with inflation. 
If we did allow the curvature to vary, our age result would 
become significantly less precise. We expect allowing w to 
vary to have little effect, as we discuss below. 

We explore the likelihood in a ten-dimensional pa- 
rameter space (six cosmological parameters plus four ex- 
perimental parameters) by Monte-Carlo generation of a 
Markov chain of parameter values as described in Chris- 
tensen et al. (2001). From the chain one can rapidly cal- 
culate marginalized one-dimensional or two-dimensional 
probability distributions for chain parameters, or derived 
parameters, with or without additional priors. Generating 
a sufficiently long chain in a reasonable amount of time re- 
quires a fast means of calculating the angular power spec- 
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trum for a given model. We describe this fast method 
briefly below and more thoroughly in Kaplinghat et al. 
(2001). 

Supernovae observations constrain the combination 
H()to better than either parameter by itself. Pcrlmutter et 
al. (1999) find for flat Universes that t = 13.0±\^(0.72/h) 
Gyr. Combining this result with our i determination 
leads to h = 0.67t;Qg, in agreement with the HST re- 
sult. Riess et al. (1998) find for arbitrary Q to t that 
to = (14.2 ± 1.7) Gyr. 

Krauss & Chaboyer (2001) estimate the age of 17 metal- 
poor globular clusters to be ice = 12.5 Gyr with a 95% 
lower bound of 10.5 Gyr and a 68% upper bound of 14.4 
Gyr. The minimum requirement for consistency, that tf = 
to — tec > is easily satisfied with a few Gyrs to spare. 
Unfortunately, the upper bound on ice is n °t sufficiently 
restrictive to set an interesting lower bound on^. 

Below we tabulate our constraints on all the model pa- 
rameters, and emphasize not only to but also With 
the inclusion of prior information on Ho, the CMB data 
provide strong evidence for f^A > 0. The same conclu- 
sion can be reached by, instead, combining the CMB data 
with observations of large-scale structure (Efstathiou et 
al. 2001) or clusters of galaxies (Dodelson & Knox 2000). 

2. METHOD 

Our first step in exploring the high-dimensional param- 
eter space is the creation of an array of parameter val- 
ues called a chain, where each element of the array, 9, 
is a location in the n-dimensional parameter space. The 
chain has the useful property once it has converged that 
P{6 e R) = N(9 e R)/N where the left-hand side is the 
posterior probability that 9 is in the region R, N is the to- 
tal number of chain elements and N(9 € R) is the number 
of chain elements with 9 in the region R. Once the chain 
is generated one can then rapidly explore one-dimensional 
or two-dimensional marginalizations in cither the original 
parameters, or in derived parameters, such as to- Calcu- 
lating the marginalized posterior distributions is simply a 
matter of histogramming the chain. 

2.1. Generating the Chain 

The chain we generate is a Monte Carlo Markov Chain 
(MCMC) produced via the Metropolis-Hastings algorithm 
described in Christensen et al. (2001). The candidate- 
generating function for an initial run was a normal dis- 
tribution for each parameter. Subsequent runs used a 
multivariate-normal distribution with cross-correlations 
between cosmological parameters equal to those of the pos- 
terior as calculated from the initial run. 

All of our results are based on MCMC runs consisting 
of 2 x 10 5 iterations. For the "burn-in" the initial 2.5 x 
10 4 samples were discarded, and the remaining set was 
thinned by accepting every 25th iteration. We used the 
CODA software (Best et al. 1995) to confirm that all chains 
passed the Referty-Lewis convergence diagnostics and the 
Heidelberger- Welch stationarity test. 

While generating the chain we always restrict our sam- 
pling to the h > 0.4 and 5.8 < z r ; < 6.3 region of parame- 
ter space. The former is a very conservative lower-bound 
on h and the latter is a simple interpretation of the spec- 
tra of quasars at very high redshift (Becker et al. 2001; 



Djorgovski et al. 2001). For some of our results we assume 
an "HST prior" which means h = 0.72 ± 0.08 (Freedman 
et al. 2001) with a normal distribution. 

2.2. Ci Calculation 

We calculate Ci rapidly with a preliminary version of 
the Davis Anisotropy Shortcut (DASh; Kaplinghat et al. 
2001). We first calculate the Fourier and Legendre- 
transformed photon temperature perturbation, A;(fc), on 
a grid over parameters u>b and Ud at fixed values of 
n k = 1 - fi to t = n A = tt* A and t = using CMB- 
fast (Seljak & Zaldarriaga 1996). From this grid, we get 
Ci for any u>b, ujd, and the primordial power spectrum 
P(k) = A(k/0.Q5Mpc^ 1 ) n by performing multi-linear in- 
terpolation on the grid of Aj(fc) and then the following 
integral: 

Ci = ^ + 2 ^° l = 8nl{l + 1) J k 2 dkAf(k)P(k). (1) 

We can get any C; in the entire model space of {ojb, ^d, t, 
Qai Slfc, P(k)} by the use of analytic relations between the 
A;(fc) for different models. For varying f^A and fife, Ci = Cf 

where I /I = s (n£, Q A )/0 s (nfc, Qa) and 9 S is the angle sub- 
tended by the sound-horizon at the last-scattering sur- 
face. For fife = 0, 9 S = s/r/o where rjo is the conformal 
time today (or, equivalently, the comoving distance to the 
horizon) and s is the comoving sound-horizon at the last- 
scattering surface. 

Altering 9 S is not the only effect of varying flh and Qa- 
Varying changes the eigenvalues of the Laplacian on 
very large scales and hence the power spectrum at the 
last-scattering surface, and both flk and £1a affect the 
late-time evolution of the gravitational potential. Both 
of these effects only affect C\ at I << 100. We therefore 
make an additional grid over the parameters ujb, k><z> f^fe 
and f^A but with smaller maximum / and k values than the 
lower-dimensional high-/ grid. The low-/ grid used for the 
calculations presented here has ranges 0.01 < Ub < 0.03, 
0.05 < uj d < 0.25, and < A < 0.85 with 4, 4, and 8 
uniformly spaced samplings of the range respectively. For 
the present application we have fixed flk = 0. For the 
high-/ grid the ranges for uib and u<i are the same but with 
twice as many samples, and A = 0.6. 

The split into a low-/ grid and a high-/ grid has been 
used by others, although for grids of Ci, not A;(fc). We 
follow Tegmark et al. (2001) in joining our grids with a 
smooth k-space kernel, g(k) = 2/(1 + exp{2k/k s ) A ) where 
k s — 1.5/s; in the integrand of Eq. 1 P(k) is replaced 
with g(k)P(k) for the low-/ grid and (1 — g(k))P(k) for 
the high-/ grid and C x = C\ ow + Cf gh . Finally, we allow 
for non-zero z ri by sending C\ lZi(z ri )Ci where lZi(z r i) 
is given by the fitting formula of Hu & White (1996). 

2.3. Likelihood Calculation 

To calculate the likelihood we use the offset log-normal 
approximation of Bond et al. (2001) which is a better ap- 
proximation to the likelihood function than a normal dis- 
tribution. We include bandpower data from Boomerang, 
the Degree Angular Scale Interferometer (DASI; Halverson 
et al. 2001), Maxima (Lee et al. 2001)) and the COsmic 
Background Explorer(COBE; Bennet et al. 1996). The 
weight matrices, band powers and window functions for 
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DASI are available in Leitch et al. (2001), Halverson ct al. 
(2001) and Pryke et al. (2001). For COBE we approximate 
the window functions as tophat bands; all other informa- 
tion is available in Bond et al. (2001) and in electronic 
form at RadPack (2001). For Boomerang and Maxima 
we approximate the window functions as top-hat bands, 
the weight matrices as diagonal and the log-normal off- 
sets, x, as zero. The Boomerang team report the uncer- 
tainty in their beam full-width at half-maximum (fwhm) as 
12. 9± 1.4 minutes of arc. We follow them in modelling the 
departure from the nominal (non-Gaussian) beam shape as 
a Gaussian. For the Boomerang Zf, Ci is therefore actually 
Ci exp(— l 2 b 2 ). Our prior for b is uniform, bounded such 
that the fwhm is always between 11.5' and 14.3'. Calibra- 
tion parameters, for example, udasi = 1 ± 0.04 are taken 
to have normal prior distributions and alter model angu- 
lar power spectra via C; — ► Up ASI C; prior to comparison 
with the reported band powers. To reduce our sensitivity 
to beam errors, we only use bands with maximum Z-values 
less than 1000. Thus we use all nine DASI bands, the first 
11 Maxima bands and all but the last Boomerang band. 

Five of the 24 DASI fields are completely within the area 
of sky analyzed by Boomerang, and three partially over- 
lap this area. We expect the resulting DASI-Boomerang 
bandpower error correlations (which we neglect) to be 
small and to have negligible effect on our results. 

3. RESULTS 

In the top panel of Figure 1 we see that neither h nor 
can be constrained well by CMB data alone. The shape 
of the contour "banana" in the top panel is determined by 
the dependence of the well-determined 6 S on h and Qa; 
lines of constant 6 S run along the ridge of high likelihood. 
Since 6 S correlates well with the age, as seen in Fig. 2, the 
ridge of high likelihood is also at nearly constant age. 

Although £Ia is poorly determined by CMB data alone 
we find that addition of the HST prior allows one to set a 
95% lower limit of Sl\ > 0.4. This same lower bound can 
be achieved by simply rejecting models with h < 0.55. 
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Fig. 1. — The posterior probability density of C2a (lowest panel) 
and contours of equal probability density in the Q\, to plane (lower 
middle panel), Qa, Os plane (upper middle panel) and Q/y, h plane 
(top panel). Contour levels are at e~ 6 - 17 / 2 , e -2 ' 3 / 2 and 0.95 of 
maximum. The lowest panel curves are for h > 0.4 (solid) and 
h = .72 ± .08 (dashed). Top panel dotted lines are at constant to- 



We can understand the to 9 S correlation with the aid 

of an approximate analytic expression given by Hu et al. 
(2001) from which we derive 

^= 0.060(2.9^ + 1.0^-1.14^) (2) 

for the fiducial values uj m = 0.15, a; a = 0.3,^6 = 0.02. 
Expanding to about our fiducial values we find 

^o = _ .i 2 f 3 .0^ + 1.2^V (3) 

Thus a change from the fiducial values by Aw m and Ao>a 
which keeps S fixed will nearly leave the age unchanged. 

In general, the parameter controlling this correlation is 
the ratio of ratios: 

(dln6 s /d\nLU m )/(dln6 s /dlnLU A ) 

(dlnio/<91nw m )/(<91n£o/dlnw A ) 
and the correlation is tightest when R = 1. R has little 
dependence on u>b- For what used to be called standard 
CDM, R = 0.75. R can be as small as 0.53 for il m = 1 and 
h = 0.72 and as large as 2.0 for n A = 0.83 and h = 0.72. 
At the maximum of the likelihood, R = 0.88(1.08) with 
(without) the HST prior. 

As a test, we have estimated the age via a direct grid- 
based evaluation of the likelihood given DASI and DMR 
data using CAMB (Lewis et al. 2000) to calculate C;'s. 
Taking the grid parameters to be to, uJd, n and A and fixing 
wdasi — 1 = — 1 = ilfe = 2 r i = wc find to = (13.6 ± 
0.6) Gyr. This agrees very well with our MCMC + DASh 
results when the same assumptions and data selection are 
made: to = (13.7 ± 0.6) Gyr. We can also reproduce the 
Netterfield et al. (2001) result, finding for Boomerang and 
DMR data (although ignoring the highest-! Boomerang 
bandpower) t = (14.32 ± 0.68) Gyr. 

In the table we show means and standard deviations 
for our ten original parameters and a number of derived 
parameters. Particularly noteworthy is 9 S , determined 
with an error of less than 3%. The agreement between 
the datasets on this number is also remarkable. From 
DASI 9 S = (0.60 ± 0.01) degrees and from Boomerang 
6 a = (0.59 ± 0.01) degrees. 
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Fig. 2. — The posterior probability density of the age of the 
Universe (lower panel) and contours of equal probability density (as 
in Figure 1) in the 9 a , to plane (upper panel). The lower panel 
curves are for h > 0.4 (solid) and h = .72 ± .08. 
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For studying the early evolution of structure, it is useful 
to know the age at redshifts in the matter-dominated era. 
For 1 << z < 100, t z x (1 + zfl 2 = Gyr = 

(16 ± 1) Gyr, where t z is the age at rcdshift z. Since MAP 1 
and Plancti 1 will determine u> m to 10% and 2% respectively 
(Eisenstein et al. 1999) they will determine t z x (1 + zf'l 2 
to 5% and 1% respectively. 

4. DISCUSSION 

Since our argument is model-dependent, it is worth 
pointing out that the model has been enormously success- 
ful on the relevant length scales (e.g. Wang et al. 2001). 
Perhaps the weakest link is the dark energy equation of 
state since we have scant guidance from observations (e.g. 
Pcrlmutter et al. 1999) and even less from theory. For- 
tunately, one can show that varying the equation of state 
away from w = — 1 at fixed 9 S has very little effect on the 
age: at 9 S = 0.6 degrees and uo m = 1/3, Ato/t = .05 Aw. 
Though we neglect the possibility of gravitational wave 
contributions to C; (see Efstathiou 2001) we do not ex- 
pect these to make much difference since 9 S is relatively 
unaffected by measurements at low t where gravitational 
waves are important. 

Our age determination has the benefit of being derived 
from observations whose statistical properties can be pre- 
dicted highly accurately using linear perturbation theory. 
We are encouraged that the observational errors are domi- 
nated by the reported statistical ones since nearly the same 
result can be derived from two independent data sets. We 
conclude that the best determination of to now comes from 
CMB data. The prospects for improving the age deter- 
mination are bright since the statistical errors (and any 
systematic ones too) will be greatly reduced by MAP data 
in the near future. 

The MCMC chains we have generated are available via 
e-mail from the authors. 

We thank M. Kaplinghat, R. Meyer and K. Ganga for 
useful conversations, B. Luey for some programming, B. 
Chaboyer for sharing results prior to publication and the 
Fcrmilab Reading Group for comments on an earlier ver- 
sion. LK is supported by NASA, NC by the NSF and CS 
by the DoE. 
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Table 1 



Parameter Bounds 



Parameter 


mean 


standard deviation 




0.021 


0.002 




0.145 


0.021 


ni 


0.49 


0.17 




6.0 


0.14 


A 


6.7 


0.56 


n 


0.96 


0.04 


MDASI 


1.00 


0.03 


^Boom 


1.07 


0.03 


^Maxima 


1.00 


0.03 


fwhm^ oom 


13'.9 


0'.3 


to (Gyr) 
h 1 


14.0 


0.48 


0.59 


0.07 




0.51 


0.17 


u>\ 


0.19 


0.11 


Mm 


0.166 


0.021 


c?jo (Gpc) 


13.95 


0.56 




35'.5 


0'.43 




168 


15 


U 


1392 


18 


H 2 


0.481 


0.024 


H z 


0.486 


0.030 




16.0 


1.0 



NOTES. — The mean and standard deviations for the ten chain parameters (top) plus derived 
parameters (bottom). For these results we use all the data with our weakest prior assumptions. 
Units of A are arbitrary. See Hu et al. (2001) for the definition of Z e q, Id, H2 and H3. The 1 
superscript indicates those parameters whose uncertainties are not well-described by a mean 
and standard deviation. For example, the posterior probability distribution for z r i is not 
significantly different from the prior one we assumed, which is uniform between 5.8 and 6.3. 



